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Abstract 

o ■ 

We study the observed enhancement of a pp system near the threshold in the process J/ip — > "fPP 
\ and e + e~ — > pp. From early studies the enhancement can be explained by final state interactions, 

which are in general taken into account with some potential models. In this work we offer a simple 
approach within quantum held theory to explain the observed enhancement. We point out that among 
different final state interactions the rescattering in a NN system though exchange of 7r is the most 
important. The effects of the rescattering is completely fixed by the well-known coupling g^NN- Our 
results show that the enhancement in J/ip — > "/pp and e + e~ — ► pp can be well described with the 
rescattering effects. 

1. Introduction 

Recently it has been reported by several experimental groups unusual enhancement of a baryon- 
antibaryon system produced near the threshold. In the decay J/ip — > ^ypp the enhancement has been 
y—i . observed by BESpQ near the threshold of the pp system. BES also has observed the enhancement in the 
decay J/ip — > K~pA near the threshold of the pA-system[2]. In the decay of B + — > K + pp and B° — > D°pp 
the enhancement has been observed by Belle near the threshold of the pp system[3]. Recently Barbar has 
\Q . reported the enhancement in the process e + e~ — > pp, AA, £°S°[U [5], respectively. These experimental 
results have stimulated many theoretical speculations^ El El HH Q21 [13]. It seems that the observed 
enhancement is a general feature for a baryon-antibayron system produced near its threshold. In this 
■ work we focus on the enhancement of pp system observed at BES and BarBar. 

The experimental data of the observed enhancement in the decay J/ip — > jpp can be described 
with an 5-wave Breit-Wigner resonance function with a peak mass at M = 1859 ± 6MeV below the 
threshold [I], and an analysis of the angular distribution of the photon suggests that the pp system is 
likely with the total angular momentum J = OpQ. Many explanations for the observed enhancement at 
BES exist. A class of explanations is that the enhancement is interpreted as the existence of a baryonium 
bound state[6j or a glueball below the threshold[7]- Another class of explanations is to take the effect 
of final state interactions into account. There are different ways to take final state interactions into 
account. One can use a complex S- wave pp scattering length |8j or use a A'-matrix formalism to include 
one pion exchange[llj. A more realistic way is by using models of NN interactions [91 [10]. These models 
are partly fixed by a well-known NN interacting potential which corresponds to the dispersive part 
of the NN interaction. The observed enhancement in e + e~ — * pp has also motivated some theoretical 
studies [lOinU [ESI US]. It is interesting to note that by taking final state interactions into account through 
models of NN interactions, the enhancement in J/ip — ► ~fpp and e + e~ — > pp can be explained |10l 115] . 
However, these models are in general complicated and contain several or more parameters which need to 
be fixed. In this work we offer a simple approach within quantum field theory to explain the observed 
enhancement in the J/tp decay and e + e~ annihilation. 

Once a NN system is produced, the final state interaction can happen in several ways. Near the 
threshold, the most important final state interaction is expected to be the rescattering chain process, the 
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rescattering can be multiple, i.e., NN — * NN ■ ■ ■ — > NN, where only NN's are in intermediate states. 
Since the momentum transfer near the threshold is small and approaches to zero, one can expect that 
the most important rescattering is through the exchange of one pion, because pion is the meson with 
the smallest mass. This can be seen from the propagator of the exchanged particle. With this argument 
one also finds that the intermediate states can only be NN states. Otherwise, the exchanged particle is 
heavier than pion and its contribution is suppressed. It is possible to have more than two particles in 
intermediate states, but the contributions from these states are suppressed by phase-space factors and 
also by propagators of particles heavier than pions. At the low energy, the coupling constant of ttNN 
is fixed as g n NN- Therefore in our approach there is only one well-known parameter. In fact such an 
effect of rescattering has been considered for the J /ip decay in [11] where the interaction is described by 
a potential. Our results are different than those in [11] and the difference will be discussed. Within our 
approach we can explain the enhancement observed in J ftp — > jpp and e + e~ — > pp. In our approach we 
will neglect the Coulomb enhancement factor, because it only has a significant effect within a few MeV 
above the threshold. 

Our work is organized as the following: in Sect. 2 we give our result for the multiple NN rescattering 
though the ttNN interaction. The effect of the multiple rescattering can be summed in amplitudes and 
analytical results for amplitudes are given. In Sec. 3 we give our numerical results for the enhancement in 
J/ip — > jpp and also for the decay of J/ip — ► TT°pp. In Sec. 4 we compare our predictions with experimental 
results of e + e~ — > pp. Sect. 5 is our summary. 



1. Rescattering in NN Systems 

As mentioned in the introduction we consider the rescattering of ./ViV system through the interaction 
between ir and ./V at low energy. The interaction is well-known and given through the effective Lagrangian 
as: 

5C = ig nNN NT l lb Nit\ N=^y (1) 

where ./V is the field for nucleon N, ir l (i = 1, 2, 3) is the pion field and T l {i = 1, 2, 3) is the Pauli matrix 
acting in the isospin space. For a simple representation of our approach and results we first study the 
rescattering process of single channel, i.e., pp — > pp- ■ ■ — > pp. Then we generalize the results of the single 
channel to those of a NN system. 

We consider a pp system produced through a vertex T, then the rescattering of the system can happen 
through the exchange of ir°, as drawing in Fig.l, where the black circle represents the vertex T. At first 
look it may be meaningless to sperate the production amplitude into a vertex part and a rescattering part, 
because the combination of the two parts is exactly the vertex for the production. However the separation 
is still meaningful if one takes the vertex in Fig.l as that obtained from a analytical continuation of the 
vertex in the space-like region into the time-like region above the threshold, where only the dispersive 
part is taken near the threshold. In that sense the absorptive part of the production amplitude comes 
only from the rescattering. We will call the vertex as "bar" vertex. 

The production amplitude in our approach is then the sum of the amplitudes in which exchanges of 
arbitrary number of 7r° happen. We denote the amplitude the rescattering though n-n exchange as T n . 
The amplitude without the exchange is denoted as To, which is given by 

T = u(pi, si)T {p \pi,p 2 )v{p 2 , s 2 ). (2) 

In the above T is the "bar" vertex. The proton carries the momentum p\ and the spin s\ and the 
antiproton carries the momentum p 2 and the spin s 2 ■ The amplitude with the rescattering through one 
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Figure 1: The Feynman diagrams for the rescattering between a proton and an antiproton through the 
exchange of ir°. The dash lines cutting diagrams are cuts. 
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7T° can be obtained with Cutkosky rules as: 

7i = 2igl NN 7T 2 J ^^S(kj - m 2 p )5{kl - m 2 p ) • ^ _ 

• ^2 u(pi, si)-y 5 u(ki, s)u(ki, s)r (p) (/ci, k 2 )v(k 2 , s')v(k 2 , s')j 5 v(p 2 , s 2 ), 

s,s' 

k 2 = Pi+P2~fa- (3) 
Similarly, one can write down the amplitude T n . The amplitude with all rescattering effects is the sum: 

oo 

% s = £ T n . (4) 

n=0 

We will work with the rest frame of the pp system. In the system the momenta are: 

I$=(E,p), p% = (E,-p). 



(5) 



We assume that the dependence of the production rate on the invariant mass q 2 with q 2 = {p\ + p 2 ) 2 
near the threshold is dominantly determined by the final state interaction or rescattering. With the 
assumption one can expand the vertex and the products of spinors in the small p: 



r (p) ( M ) 

v(k 2 ,s)j 5 v(p 2 ,s 2 ) 



T^(E,p) = T^\m p ,0){l + O(P)} 

V Hs)*-(k 2 -p 2 )ri(s 2 ){l + 0({3)}, 
2 



Am p 



q 2 = (pi +P2) 2 , 



(6) 



where ^ is the two-component spinor for the proton and 77 is that for the anti-proton. The above expansion 
implies that we take the nonrelativistic limit (3 — > 0. In the limit we can classify the state of system with 
the sum of spins s = s± + s 2 by using: 
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£,(siW(s 2 ){si, s 2 \s = 0,s z = 0) = 



^(si)r7 t (s 2 )(si,S2|s = Mz) 



(7) 
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where / is a 2 x 2 unit matrix, (si, S2\S, S z ) is the Clebsch-Gordan coefficient and e(s z ) is the polarization 
vector for the case s = 1. We can decompose the amplitude into a s = 0-part and a s = 1 part: 

Ti = 7i a=0 + Ti i=1 , 



= i9l NN * 2 j ^(kl-ml)5(ki-ml) 



•Tr [<r • (ki - pi)<7 ■ (pi - ki)] A p p, 



7T 1 = i 9 l NN ^ 2 j ^^l-rnl)5{kl-rnl) 



•Tr [<r • e*(s z )cr • (ki - pi)er • e(s'Jo- ■ (pi - ki)] B pP (s f z ), 
A pP = {0,0\s 1 ,S2)u(po,s 1 )r ( - p \mp,0)v(po,s 2 ), 
Bpp{s z ) = (l,s z \si,s 2 )u(po,s 1 )T ( - p \mp,0)v(p ,s 2 ). (8) 

The summation over repeated spin indices is implied. The amplitude A p p is the amplitude for the pp in 
the state 1 S'o before the rescattering, the amplitude B p p(s z ) is that for the pp in the state 3 Si before the 
rescattering. 

It is straightforward to obtain the result of T\ by performing the phase-space integral. We have for 
the s = part: 

Tr° = % 2 ) = %^(2- y ln(l + ^)), y =M. (9 ) 

With arguments of symmetry one can show that the total spin s will not be changed after multiple 
scattering. Hence we have: 

rr- q~S=0 _|_ SJ-S = 1 q~S=0 \ 1 q~S=0 q~S = l \ ^ q~S = l / -1 

J-rs — I rs ~r 1 rs i 1 rs ~ / , J n i 1 rs ~ / t 1 n • V iu 7 

n n 

Inspecting the structure of the amplitude T s=0 one easily finds the result and the sum: 



oo 



T 



s=0 



idtf) ^ App, r r r = E^r = j-f^y. (ii) 



From the above expressions one can see that the g 2 -dependence in the amplitude appears through the 
variable y combined with in Eq.(9). Hence the energy scale m w characterizes the g 2 -dependence. In 
the limit (3 — > the amplitude is proportional to m~ 2 . If other particles are exchanged, then instead of 

their masses appear in y and characterize the g 2 -dependence. Since these particles must have masses 
M larger than m^, their contributions to the amplitude will be suppressed by m 2 /M 2 in comparison 
with that of 7r-exchange and lead to a small correction to the g 2 -dependence in the above. This is the 
argument given in the introduction to support our approach. 

Now we turn to the s = 1 part. In our approach, the pp system produced from the vertex without 
scattering has the orbital angular momentum I = 0. After the rescattering with the interaction in Eq.(l), 
the orbital angular momentum can be I = and 1 = 2. The total angular momentum J remains the 
same as J = 1. With the rescattering the amplitude will have two components, one is with i = 0, another 
is with 1 = 2. The two components are at the same order of (3. We introduce the notation for the two 
components: 

{pp)o(s z ,m z ) = e*(s z ) ■ e(m z ), (pp) 2 (s 2 ,m 2 ) = p • e*(s 2 )p • e(m z ) - ^ e *( s z) ■ e(m z ). (12) 
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With the notation we can present the result for 7| as: 



T/ =1 =B pP {m z ) (1,0) 
where M is a 2 x 2 matrix given by: 
M = 



64vr 



(pp)o(s 2 ,m 2 ) 
{pp)2{sz,m z ) 



(13) 



-Moo 






-M 2 2 


(- 




")■ 







-§(2-yL w ), (2 - 6y + (3y 2 + 2y)L y ) 

(l-2y + (y* + %)L y ), (-§ - 2, + (y 2 + f ) L y ) 



(14) 



In the above the row vector (1,0) indicates that the pp pair produced by the vertex is with £ = 0. 
The matrix elements M.20 and M.22 are actually identified through the calculation of T2. They will not 
contribute to T\ here. The above result also indicates that the rescattering of the pp pair with s = 1 is 
already a coupled channel problem. Again, the structure of 7^ =1 can be easily found and the sum can 
be obtained: 



Tn' 1 = B pP (m z ) (1,0) 



NN 



64vr 



M 



T 



s=l 



Y t Tr 1 = B vP {m z ){l,Q) 

n=0 



(pp)o(s z ,m z ) 
(pPh(s z ,m z ) 
~i -1 



M 



NN 



64vr 



M 



(pp)o(s z ,m z ) 
(pph(s z ,m z ) 



(15) 



The above results can be generalized to the rescattering of a NN system with the given interaction by 
taking the isospin factor into account. We will neglect the mass difference between protons and neutrons 
and that between different 7r's. These differences are small and only give small corrections to our results. 
For T{ the isospin factor can be determined by considering the scattering NiN^ — > NjNi through one-7r 
exchange. The isospin factor in the amplitude is: 

{T a )ji{T a ) k i = 28^8* - SjiSu. (16) 

We denote the "bar" vertex for pp and nn as and T^, respectively. We have for the s = part: 



T 

-'n 



s=0 



T 

rs 



s=0 



(*A pp ,*A nn ) id(q )C 

00 

%i = {Aw 



n=0 



{PP) ) 

(nn) J ' 
I - id(q 2 )C 



(pp) 
(nn) 



with 



C 




(17) 



(18) 



C is a 2 x 2 matrix acting in the isospin space, (pp) and (nn) denotes the 1 S'o states of pp and nn, 
respectively. ^l nS is obtained by replacing in A pp with r( ra ). We note for the above result that the 
matrix C is diagonal if we use the basis of isospin. 
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For the s = 1 part the results are: 



J-S 
x 1-1 



rs 



(B p p(m z ),0,B nn (m z ),0) 



+i 



647T 



£ IT 1 = (Bpp{m z ),0,B n n(m z ),0) 

n=0 



E-i 



( (pp)o(s z ,m z ) \ 

(pp)2(s z ,m z ) 

(nfi) (s z ,m z ) 
\ (nn) 2 (s z ,m z ) J 

I (pp)o(s z ,rn z ) 
(pPh(s z ,m z ) 



64vr 



C tg> M 



-i 



nn) (s z ,m z 
\ (nn) 2 (s z ,m z ) ) 



(19) 



where E is a 4 x 4 unit matrix. Now we can use the above results to make predictions for relevant 
experiments. 

Before turning to our numerical predictions, it is important to know the g-region where our approach 
is applicable. In our approach we have considered the rescattering effect among a NN system which 
contains the two most lightest baryons. The next lightest baryon is A. In the range of q 2 > 4m^ A's 
will enter into the rescattering as intermediate states which is not included in our results. Therefore our 
result may not be useful for q 2 > Am\ = (2 • 1115MeV) 2 . We have also taken the nonrelativistic limit 
(3 — > in our approach. The effect of the next-to-leading order in (3 can be important if (3 is large. For 



2m p ~ 



200MeV (3 takes the value ~ 0.2 which is not large. Therefore we expect that our results 
should be applicable in the range ^/q 2 — 2m p < 100 ~ 200MeV. It should be noted that it is possible to 
include higher-order effects in (3 and other neglected effects discussed in the above. 



3. The Enhancement in the Decay J /ip — > jpp 

There is an enhancement of the pp system near threshold in the decay J/ip — > jpp observed by BES 
and the experimental results favor that the pp pair near the threshold is in the state 1 Sq. This corresponds 
to our result with s = 0. In the experimental result of BEST], most data points lie in the region with 
q 2 > 4m 2 . Therefore nn pairs can appear in the intermediate states. It is expected that not only a pp 
but also a nn can be produced through the "bar" vertex in the radiative decay of J/tp. From Eq. (17) 
we obtain the decay amplitude with the rescattering effect for q 2 > &m 2 n : 

rys=o -App + ■A.nn 1 ~\~ id(q 2 ) App ■A-nn 1 3id(q 2 ) /on} 

rs " 2 {l-id{q 2 )) 2 +Ad 2 {q 2 ) 2 (1 - id{q 2 )) 2 + U 2 {q 2 ) ' 1 ' 

where the first term is with the isospin 1 = and the second is with 1 = 1. The constants A p p tn n are 
unknown because the isospin in the decay is violated in general because of the electromagnetic interaction. 
In our numerical predictions we will take the constants A p p t nn as free parameters to fit the BES results. 

In Fig. 2 we draw the BES date from [1\ and our fitting results. We have fitted the experimental date 
in the region with \[q 1 ~ — 2m v < 150MeV. The fitted parameters are A p p ~ 0.895 and A n n ~ 1-027 with 
X 2 /d.o.f. ~ 1.14. From Fig. 2 we can see that our predicted q 2 -dependence in Eq.(20) matches fairly well 
with experimental data in the region with 150MeV> ^/q 2 — 2m p . Our prediction fails to describe the 
experimental data in the region with \f(p — 2m p > 150MeV. The reason can be that the higher-order 
effects in (3 or other neglected effects can become important in this case. However, from Fig. 2 we see that 
the observed enhancement happens in the region with \f~q 1 ~ — 2m v < 150MeV, we can conclude that the 
enhancement is fairly well in agreement with our result. Also, our results indicate that the pp near the 
threshold is almost in the state with 1 = because A n n + A p p » A n n — A p p from our fitting results. 
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Figure 2: Comparison of BES results J/tj) — > jpp with our predictions. 

In [TT] the rescattering effect in the 1 = channel through 7r-exchange are studied where the interac- 
tion is given by a potential. We note that the potential used in |Tl] is different than the potential derived 
from the interaction Lagrangian 5C in Eq.(l). This leads to the difference between our predictions and 
those in [11] . Numerically, the results for the spectrum obtained in [11] look very similar to ours in the 
1 = 1 channel, which we also give in Fig. 2. From Fig. 2. it is clearly that the enhancement can not be 
explained only with the rescattering effect in the 1 = 1 channel. Our result is also different than that in 
[9], in which the enhancement is explained with the pp system in the 1=1 state. 

In experiment the decay J/ip — > ir°pp has been studied, but no enhancement near the threshold of the 
pp system has been observedpQ. It has been a puzzle. In general the final state interaction or rescattering 
is more complicated than the radiative decay, because one has three hadrons in the final state and the 
final state interaction exists not only in the pp system but also between the system and 7r°. If we assume 
that the rescattering effect of the NN system is dominant, then the observation can also be explained 
in our approach. In the decay the pp system is in the 1 = 1 state. Near the threshold the pp system 
can only be in a 3 5i state or a 1 Pi state. For the P-wave state the decay amplitude is suppressed by 
an extra j3 in comparison with that of a 5-wave state. If we take the pp system produced through the 
"bar" vertex as a 3 Si state, we can give predictions in our approach. In Fig. 3 we give our numerical 
results from Eq.(19) in comparison with BES data. It can be seen that our approach can explain the 
(^-dependence fairly well near the threshold. 

Recently, BES has reported that no enhancement in the decay J/tp — > ujpp has been found near the 
threshold [17]. In this decay the pp system is in the state with 1 = 0. If one assumes that the rescattering 
effect of the NN system is dominant as in the above, one expects from our approach that there will be 
an enhancement. The existence of the enhancement is also predicted in [18]. However, the assumption 
may not be correct. In general, for a complete prediction in this case one needs to work out rescattering 
effects among three hadrons. This deserves a further study. 
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Figure 3: Comparison of BES results of J/ip — * ir°pp with our predictions. 



4. The Enhancement in e + e — > pp 

The total cross section of e + e~ — > 7* - 

a(e + e~ — ► pp) 



pp can be expressed in terms of form factors of protons: 
4vra 2 /3 



3q 2 



G 



(P)|2 
M I 



+ 



2m P |^(p)|2 
n 2 \ U E I 



21) 



where Gm and Ge are the Sachs form factors of protons. They are related with the Dirac- and Pauli 
form factor. The definitions of these form factors are: 



(P(P1)P(P2)|^|0) 

Gm 



u(pi 



.0 



/Ml 



^F^+C—q^q 2 ) 



Fi + F 2 , G E = F 1 + 



4M 2 



q = Pi +P2- 



(22) 



In our approach the measured form factors are related to the "bar" form factors through rescattering. 
Since we are interested in the behavior near the threshold, we expand the above spinors in (3. It should 
be noted that with our formula given in Eq. (19) it is possible to keep the relativistic effect from the 
spinors. Including this effect we have for the matrix element: 



(p(pi)p(p 2 )|J|0) 



-2M (F^q 2 ) + F 2 (q 2 



£*<TTf 



2M 2? 



1 

a pp - -a 



M 



~~ 1 ' 

(7 ■ pp (T 

^ 3 



IdI 2 

V-m(q 2 )^ ] crr 1 + - 



(23) 



where • • • stand for contributions at /? 4 and are neglected in our numerical predictions. We denote the 
form factors in the "bar" vertex, i.e., the "bar" form factors, as F\ 2 and the inverse of the 4x4 matrix 



S 



in Eq.(19) as: 



A 



Pal 



NN 



32(2tt) 2 

then we have the form factors of protons with the rescattering effect near the threshold as: 



(F 1 (q 2 ) + F 2 (q 2 )) 



+ F 2 {q 



2, IPl 



3M 2 



IPl 



2M 2 



F 1 (q 2 )-F 2 (q 2 



(p) 



(p) 



An (Fi + F 2 ) iP) + A 31 (Fi + F 2 ) {n) 
An (Pi + F 2 ) (P) + A 32 (Fi + F 2 ) {n) 



(24) 



(25) 



where the last equation comes from the 1 = 2 part. The indices (p) and (n) denote quantities for proton 
and neutron, respectively. In the above the "bar" form factors should be taken with q 2 = Am 2 in our 
approach, see Eq.(6). Therefore, the g 2 -dependence of the form factors and of the cross-section comes 
from the rescattering effect and the relativistic correction in the expansion of spinors. It should be noted 
that there will be only one form factor if the relativistic effect is neglected, i.e., Ge = Gm- 
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Figure 4: Comparison of Barbar results with our predictions. 

With the above results one can obtain predictions for form factors and the cross-section near the 
threshold. If one knows the "bar" form factors, the form factors are completely determined. As discussed 
in Sect. 2., these "bar" form factors should be obtained as an analytical continuation of form factors 
from the space-like region with q 2 < to the time-like region with q 2 = Am 2 . There are intensive 
studies of experiments and in theory for the form factors in space-like regions. However, the continuation 
with existing results may be problematic, since in the continuation q 2 will cross the unphysical region 

< q 2 < Ami. We wil1 take the combination g£ } = (Pi + F^ and G$ = (f x + F 2 ^ P) as two free 
parameters to fit the Barbar results. To see the enhancement clearly, one can use the effective form factor 
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which is measured in [3]: 

.(e+E-^p) n= +gZ( 1+ *!£\. (26) 

In the Fig3. we give our numerical results where the Barbar results in [3] are also plotted. Our results 
are obtained by fitting the BarBar data in the region with q 2 < (2.230GeV) 2 with the fitting results: 

0.213, 0.427, X 2 /d-o.f. « 1.02. (27) 

From Fig. 4 we can see that the behavior of the enhancement is well described with our results. From our 
result given in Fig. 3 for J/tp —* ir°pp, where the pp is with 1=1, one may expect that the enhancement 
in e + e~ — ► pp is mainly due the enhancement in the channel with 1 = 0. But from our results we can see 
that a substantial 1 = 1 component is needed to match the experimental results. 



5. Summary 

In this work we have studied the enhancement of a pp system near the threshold in the process J ftp — > 
^ypp and e + e~ — ► pp. From early studies the enhancement can be explained by final state interactions, 
which are in general taken into account with some potential models. We have given arguments to support 
that among these interactions the rescattering effect of NN system though exchange of tt is the most 
important. Hence we have proposed a simple approach in the framework of quantum field theory with 
the well-known ttNN interaction to explain the enhancement, where the rescattering effect is completely 
fixed by the well-known coupling q^nn- It turns out that the enhancement in J/ip — * jpp and e + e~ — > pp 
can be well described within our approach. In these two cases the final state strong interaction can only 
happen between the NN system. With the assumption that only the rescattering effect in the NN 
system is dominant in J/ip — ► ir°pp, we can also explain why there is no enhancement in the decay. 

We have only studied in detail the enhancement with the final state which contains a pp pair as 
hadrons. For other cases like the enhancement in 5-decay and the annihilation of e + e~ into baryon pairs 
other than a pp pair, extensions of our approach are needed. This will be studied in a future work. 
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